Facet ridge end points in crystal shapes 
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Equilibrium crystal shapes (ECS) near facet ridge end points (FRE) are generically complex. 
We study the body-centered solid-on-solid model on a square lattice with an enhanced uniaxial 
interaction range to test the stability of the so-called stochastic FRE point where the model maps 
exactly onto one dimensional Kardar-Parisi-Zhang type growth and the local ECS is simple. The 
latter is unstable. The generic ECS contains first-order ridges extending into the rounded part of the 
ECS, where two rough orientations coexist and first-order faceted to rough boundaries terminating 
in Pokrovsky-Talapov type end points. 

PACS number(s): 64.60 Fr, 68.35 Bs, 64.60 Ht, 68.35 Rh 
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Equilibrium crystal shapes (ECS) have been studied 
over the span of many years jl| . Many features of their 
thermal evolution are theoretically well understood and 
have been observed experimentally. For example, sur- 
face roughening of facets are realizations the Kosterlitz- 
Thouless transition ||l| . They have been observed in He- 
lium crystals 0, metal surfaces Q, ionic solids [Q and 
organic crystals ||5|]. Smooth boundaries between flat 
facets and rounded rough regions, like in Fig. |l|(a), be- 
long to the so-called Pokrovsky-Talapov (FT) universal- 
ity class. They have been observed in, e.g., lead [g] and 
Helium crystals [[7|. An aspect for which experimental 
data is still sparse is the ECS structure near facet ridge 
end point (FRE). For example, in NaCl crystals, facet 
ridges seem to vanish very quickly as temperature is in- 
creased [B|. Many other experiments reproduce growth 
shapes, where FRE points are not represented properly 
(rough surface orientations are suppressed because they 
grow faster than facets). 

Lately the ECS structure near FRE points has gained 
theoretical attention, with the discovery of a link 
with Kardar-Parisi-Zhang (KPZ) type non-equilibrium 
growth Q. Fig. [|(a) represents the simplest and un- 
til now canonical ECS structure near FRE points. The 
sharp ridge between two faceted orientations terminates 
and splits into two PT facet-to-round boundaries. This 
structure is realized in the so-called body-centered solid- 
on-solid (BCSOS) model with next-nearest neighbor in- 
teractions only pO| , which is equivalent to the exactly 
solved six-vertex model. The FRE point where the PT 
lines join the facet ridge is very intriguing. At this point, 
the model can be mapped exactly onto a one dimensional 
BCSOS growth model that lies in the KPZ universality 
class 0]. To be more precise, the transfer matrix of the 
two dimensional equilibrium model maps onto the mas- 
ter equation of the one dimensional growth model. The 
spatial direction parallel to the facet ridge plays the role 
of time in the dynamic process. 

This seems to settle the issue, but the FRE point in this 
exactly soluble model has too much symmetry (by being 
stochastic) to leave us confident that its scaling behav- 



ior is universal and that the ECS of Fig. H(a) represents 
the generic shape. Recently, we performed a study of 
the ECS in a BCSOS model on a honeycomb lattice pi[ . 
There we did not find the structure of Fig. |l|(a), but in- 
stead shapes equivalent to Fig. |l| (b) and (c). In (b) the 
facet ridge extends into the rough region as a first-order 
rough-to-rough (FOR) line. It terminates at a critical end 
point. This possibility had been suggested before iQ but 
had not been observed in any model beyond mean field 
theory. Our result was somewhat unsatisfactory however 
since the FOR line remained extremely short. In another 
part of the phase diagram we found an ECS similar to 
Fig. |l|(c). The boundaries between the faceted and the 
rounded rough regions are partially first-order. These 
segments became quite long and we could establish the 
scaling behaviour of the critical end points (PTE points) 
where the transition becomes PT like. PTE scaling is in 
agreement with a simple theory where the free energy is 
assumed analytic on the rough side and expandable in 
powers of the tilt angle [0-0[. PTE points have been 
observed in Si(113) flj], but with scaling exponents dif- 
ferent from the above mean field type values. This can be 
attributed to long range attractive step-step interactions 
in the context of step bunching |l^,|l^ . 

In this letter, we present a study of the BCSOS model 
on the square lattice with enhanced interaction range. 
We address directly the stability of the stochastic FRE 
point. The FOR lines are well developed in this model. 
They appear also as stand-alone ridges inside the rounded 
parts of the crystal, as shown in Fig. |l|(d), representing 
spontaneous symmetry breaking within the rough phase 
(coexistence between two rough surface orientations). 

Consider a square lattice, oriented diagonally. The 
most general Hamiltonian with only next-nearest neigh- 
bor interactions, the one that is exactly soluble, can be 
written as 
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with X + 2,y {x,y + 2) and x,y next nearest neigh- 



bors in the horizontal (vertical) direction and hx,y = 
0;±1,±2,--- the height variables. Nearest neighbor 
heights must differ by one. The fields Eh and Ey couple 
to the slopes, which we denote as Qh and Qv respec- 
tively. In this diagonal set up, the time like direction 
of the transfer matrix coincides with the facet ridges in 
Fig. ^ The ECS of this model takes the structure of 
Fig. [l|(a) when exp(— Lj:/) > exp(— Ly) + 1, in particular 
when the step-step interactions are attractive in the hor- 
izontal direction and repulsive in the vertical direction. 
The FRE point has the stochastic KPZ character men- 
tioned above ||l^,^ . To study the stability of the stochas- 
tic FRE point we need to break the special symmetries 
of Hq. Increasing the interaction range will do this. We 
choose uniaxial interactions for numerical convenience. 
The simplest is to add interactions between next-to-next 
nearest neighbors lying in the same horizontal row 



H = Ho + Y,i^l2 5{K^y, K+i,y ± 2) 

x.y 

+Mi S{hx,y, hx+4,y ± 4)] 
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with S{x) Kronecker delta's. This involves two coupling 
constants, M2 and M4, because those further neighbors 
can differ in height by 0, ±2, or ±4. The details of the 
numerical analysis are the same as in our previous pa- 
per [IJ] and outlined there in detail. The eigenvalues of 
the transfer matrix yield the exact free energy as func- 
tion of tilt angle, /(EvtQhjN), at all temperatures in 
a semi-infinite strip geometry of width N. The ECS fol- 
lows from this by a iV — > 00 finite size scaling analysis, 
combined with a Legendre transform. The latter is in 
essence the Wulff construction Q . 

The phase diagram at Lh = —0.69 and Ly — 0.69 for 
zero tilt fields. Eh = Ey = 0, is shown in Fig. g. This 
is a representative example. In the upper right corner, 
where M2 and M4 are both large and positive, the sur- 
face is flat. As M4 is decreased, the surface undergoes 
a first-order transition to a faceted phase. The slope of 
the surface changes abruptly from (0, 0) to two coexisting 
{Qh,Qv) = (iliO) facets. As A/2 is decreased, starting 
from the flat phase, the surface first roughens, via a KT 
transition and then enters a spontaneously tilted rough 
phase U%. This is a rough-to-rough faceting transition 
with yet unknown anisotropic scaling properties. In the 
tilted phase, two rough phases coexist. The rounded 
part of the ECS contains a first-order ridge as shown in 
Fig. 0(d). The jump in orientation is caused by the com- 
petition between the L and M interactions, and changes 
continuously throughout this phase. It locks in smoothly 
to Qh = ±1 at the FT boundary into the faceted phase. 

The dashed line in Fig. inside the faceted phase does 
not represent a phase transition. It marks the boundary 
between the two different shapes, Figs. ^ (b) and (c). In 
region (c), the edge between each facet and the rounded 
rough region is sharp near the FRE point. The scaling at 



the PTE point (where the edge becomes a conventional 
PT boundary) is consistent with our previous study, i.e., 
with an analytic free energy functional in terms of the 
surface tilt angle. The correlation length in the direction 
along the facet ridge scales compared to the perpendicu- 
lar one as ^|| ~ l^ with z = 3 O-O. 

In region (b) the first-order facet ridge boundary be- 
tween the (±1,0) facets extends into the rounded re- 
gion as a FOR line, as shown in Fig. iVb). Unlike 
our previous study on the honeycomb lattice |11[ , the 
FOR line becomes quite long and is unambiguously re- 
solved. For example, at point (L/f , _Ly,M2,M4) — 
(-0.69,0.69, -0.02,0.02) the end point of the FOR line 
lies at Ev = 1.04 ± 0.01, while the location of the FRE 
point is known exactly, Ey = 0.908. At the phase bound- 
ary into the spontaneously tilted rough phase, the ECS 
changes from Fig. ^c) to Fig. [|(d). In the latter the 
facet-to-round boundary never touches the FOR line. 

The scaling behaviour at the FOR end point is an im- 
portant issue. A naive guess is to assume the free en- 
ergy is analytic in terms of the tilt angles. It should 
then take the asymptotic form f{QH,Ey) = A{Ey — 
Ey)Q^ + BQjj, such that at the FOR critical end point 
the free energy scales like / ~ Qj^. The free energy gap 
between the Qh = and Qh = ^/L sectors would be- 
have as L[f{l/L) — f{0)] ^ L^^, which translates into an 
anisotropic scaling exponent z = 3. The latter is inconsis- 
tent with our numerical data. Fig. a illustrates the finite 
size scaling estimates for z at two FOR end points, lo- 
cated respectively at Ey = 1.05 for {Lh, Ly, M2, M4) — 
(-0.69,0.69,-0.02,0.02), and at Ey = 0.62 for 
{LH,Ly,M2,M4) = (-0.41,0.69,0.05,-0.05). We mea- 
sure the largest and second largest eigenvalues, Aq and 
Ai, in the transfer matrix sector with periodic boundary 
conditions. The gap m = ln(Ao/Ai) scales with system 
size like m ~ N~^. Fig. |3| illustrates that the finite size 
scaling values of z are rather sensitive to the estimate 
of the location of the endpoint, and that corrections to 
scaling vary strongly throughout the phase diagram, but 
z is clearly smaller than 2. Therefore, the above naive 
description is incorrect. 

A second natural guess for the scaling properties of the 
end points of FOR lines is that they belong to the KPZ 
universality class. The transfer matrix is not stochastic, 
but maybe some of its properties survive. The numerical 
data in Fig. ^ do not exclude the KPZ value z = 1.5. 
To test this possibility in more detail, we check for a dis- 
continuity in surface orientation, AQy, at the end point 
in the direction along the FOR line. Qy jumps at the 
stochastic FRE point in Fig. H(a). This represents the 
growth velocity in 1 + 1 dimensional KPZ dynamics. At 
the end points of the FOR lines, on the other hand, this 
discontinuity vanishes as function of system size N, like 
AQy r^ N~'^ with X ~ 2.2 ±0.1 at all points we checked. 
The discontinuity in surface tilt along the FOR line in 
the perpendicular direction, AQh, acts as order param- 



eter. It vanishes continuously, probably as a powerlaw, 
AQff ^ \Ev — Ey\^. The discreteness in strip width N 
allows only commensurate rational values for Qh- This 
makes an accurate scaling analysis impossible. However, 
we are confident that 0.3 < /3 < 0.5. 

In Figs. |l] (b) and (c) the FRE point has first-order 
characteristics. The crossover region in between these 
structures, marked by the dashed line in Fig. 0, is the 
only place we can hope to find Fig. Wa) type structures. 
However, this crossover region is actually blurred and 
contains more complex crystal shapes. At the stochastic 
points M2=A/4=0 the crossover is simple. On approach 
from the (b) side, the FOR line shrinks in length and 
vanishes completely, followed by the emergence of PTE 
points out of the FRE point on the (c) side. 

On approach of the crossover region below the stochas- 
tic point, marked 1 in Fig. g, the FOR line shrinks 
rapidly as well, but just before it disappears completely, 
the PTE points start to emerge already along the two 
facet-to-round boundaries. The intermediate ECS shape 
is a superposition of structures (b) and (c) , as shown in 
Fig. ^(a). Fig. || shows the free energy as a function of 
surface tilt Qh at Ey = 0.69 for (L^, Ly, M2, M4) = 
(0,0.69,-0.14,-0.40). The crucial feature is the simul- 
taneous presence of a central maximum at Qh = fol- 
lowed by two nearby minima (this creates the FOR line) 
and concavity of the curve at the Qh = ±1 edges (the 
PTE points). We refer to our other paper for a detailed 
discussion on how such features translate into ECS prop- 
erties |ll|. The choppiness in Fig. is due to the limi- 
tation to rational values of Qh (multiples of 1/N). The 
numbers are known to arbitrary accuracy and their finite 
size scaling converges smoothly. 

This superposition of the FOR and PTE shapes is pre- 
dictable, because it is hard to imagine how local fea- 
tures oi /{Qh) at opposite tilt angles can be topologically 
linked to each other, except by special symmetries, like 
stochasticity The crossover region is very narrow, but 
numerically significant. For example, in Fig. [s] the FRE 
point lies at {Eh,Ev) = (0,0.668 ± 0.004) and the FOR 
end points at (0, 0.70 ± 0.01), compared to the virtually 
exact location (±0.15808,0.87105) of the PTE points. 

Above the stochastic point, marked 2 in Fig. ||, the 
crossover is more intricate. Here, the free energy as a 
function of Qh develops extra wiggles and a slightly con- 
cave part at intermediate Qh- On approach of the dashed 
line, the FOR line shortens dramatically but then splits, 
see Fig. 0(b). The two FOR end points rapidly move 
outwards to merge with the facet-to-round boundary into 
PTE points. 

At the stochastic point the free energy curve /{Qh) is 
completely featureless. It is flat. All tilt angles are de- 
generate, due to the special symmetries associated with 
stochasticity. In region 2 of the crossover region /{Qh) 
becomes almost flat, but retains some wigglyness. Some 
of this is thermodynamically unstable (skipped by the 



Legendre transform), but the ECS is susceptible to com- 
plex hairy behaviour, instead of simplifying into the fea- 
tureless structure of Fig. ||(a) . Complexity is indeed what 
we find, although the splitting remains ambiguous, be- 
cause the splitted FOR lines are extremely short. We 
can not rule out that their stability is a finite size scal- 
ing artifact. Therefore, we determined also the (effec- 
tive) anisotropic exponent z, in region 2 pretending the 
ECS takes the simple form of Fig. 0(a). The results are 
consistent with the KPZ value z — 1.5, but this can be 
attributed to crossover scaling from the stochastic point. 

In conclusion, the previous canonical picture of equi- 
librium crystal shapes near facet ridge end points where 
the facet ridge splits into two PT lines, like in Fig. |^(a), is 
too naive. We established that the stochastic FRE point 
is unstable and for the first time observed FOR lines un- 
ambiguously. Generic experimental equilibrium crystal 
structures near FRE points will contain FOR lines stick- 
ing into the rough rounded parts and first-order facet-to- 
round segments, Figs. [^ (b) and (c) and Fig ^. We find 
also stand-alone FOR lines, ridges inside the rounded 
part of the crystal where the surface orientation jumps, 
as in Fig. |l|(d). The scaling properties of the critical 
end points of those ridges, and the critical point at their 
birth (the rough to tilted phase boundary in Fig. ||) need 
further study. This research is supported by NSF grant 
DMR-9700430. 
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FIG. 1. Equihbrium crystal shapes in the BCSOS model 
with enhanced interaction range: (a) ECS in the exactly sol- 
uble square lattice BCSOS model with stochastic FRE point. 

(b) ECS with a first-order line extending into the rough area. 

(c) ECS with first-order facet-to-round boundaries and PTE 
points, (d) ECS with a spontaneous tilted rough phase, i.e., 
with a first-order ridge inside the rough phase. 
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FIG. 2. Phase diagram of the model of Eq. |2| for 
Lh = -0.69, Lv = 0.69. The labels (b), (c), and (d) re- 
fer to the ECS in Fig. |i|, and 1 and 2 to Fig. |(a) and (b). 




FIG. 3. The effective dynamic exponent as a function 
of strip width A'' at the end points of two FOR line: 
(a) {Lh,Lv,M2,M4) = (-0.69,0.69,-0.02,0.02) and (b) 
{Lh,Lv,M2,M4,) = (-0.41,0.69,0.05,-0.05). The dashed 
lines represent the uncertainty in the location of the end 
points. 
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FIG. 4. The ECS near the FRE point in the crossover re- 
gion between the shapes of Fig. |l| (b) and (c). 
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FIG. 5. The free energy as a function of Qh for 
(Lh, Lv,M2, Mi) = (0.0, 0.69, -0.14, -0.40), at Ev = 0.683. 
Only the Qh > side is shown. fiQn) is mirror symmetric. 



